Quantum Nucleation in a Ferromagnetic Film Placed in a 
Magnetic Field at an Arbitrary Angle 



Gwang-Hee Kim* 
Department of Physics, Sejong University, Seoul 143-747 , Korea 
(Received ) 

Abstract 

We study the quantum nucleation in a thin ferromagnetic film placed in a 
magnetic field at an arbitrary angle. The dependence of the quantum nucle- 
ation and the temperature of the crossover from thermal to quantum regime 
on the direction and the strength of the applied field are presented. It is found 
that the maximal value of the rate and that of the crossover temperature are 
obtained at a some angle with the magnetic field, not in the direction of the 
applied field opposite to the initial easy axis. 
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Consider a ferromagnetic film of thickness Zt and its plane is perpendicular to the easy 
axis determined by the magnetocrystalline anisotropy energy depending on the crystal sym- 
metry. Since the magnetic anisotropy energy has a time reversal symmetry in the absence 
of an external magnetic field, M has at least two easy directions. By applying a magnetic 
field H in a direction between perpendicular and opposite to the initial easy axis along 
which the system is magnetized, the symmetry of two equivalent orientations of M is bro- 
ken. According to the classical theory of nucleation, the switching of the magnetization 
occurs for a magnetic field greater than the critical field He at which the spins of a saturated 
ferromagnet cease to be aligned. Besides the magnetic anisotropy energy Ea-, the exchange 
energy E^x and the magnetic dipole-dipole interation E^d determine the switching process 
of the system. In this paper we simplify the system with small E^d compared with Ea and 
Eex- At finite temperature, even for the field less than there is a nucleation due to 
thermal fluctuation whose rate is proportional to exp{—U/kBT) where U is an energy bar- 
rier related to the magnetic anisotropy energy, the external magnetic field and the exchange 
energy. At a temperature low enough to neglect the thermal activation process, the mag- 
netic nucleation which is independent of temperature, may occur due to quantum tunneling. 
For the dynamical situation, it is important to include the effect of the environment on the 
quantum tunneling rate caused by phonons, 0] nuclear spins, p and Stoner excitations and 
eddy currents in metallic magnets. However, many studies have shown that while these 
coupling might be crucial in macroscopic quantum coherence, they are not strong enough 
to make the quantum tunneling unobservable. 

The question of the quantum nucleation were studied by Privorotskii |^ who suggested 
the WKB exponent based on the dimensional analysis. Later, Chudnovsky and Gunther 
proposed the quantum nucleation in a ferromagnetic thin film at zero temperature by 
using the instanton technique. In their work it was considered the system with the external 
magnetic field applied in the opposite direction to the initial easy axis, and presented the rate 
of quantum nucleation which allows experimental observation with a sufficient probability. In 
this paper we shall extend their considerations to the system with the magneic field applied 
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in a range of angles 7r/2 < < vr where 9h is the angle between the initial easy axis (z-axis) 
and the external magnetic field, which, in its method and spirit, is the generalization of Ref. 
P] to problems with a non-uniform rotation of M. 

Using the parameter defined to be e = 1 — H/Hc, at zero temperature the nucleation 
exponent Bq oc e^/^ in a wide range of angles, not e^/^ which was obtained by Chudnovsky 
and Gunther for 9h = t^- [@] The crossover temperature for quantum nucleation to ther- 
mal activation has a weak dependence on e, Tc oc e^^'^ in a wide range of angle [§, compared 
with e^/'^ at 9h = tt. [§] Also, it is found that the direction of H which favors the quantum 
nucleation lies in the range of angle /2 < 9h < 'n: ■, not tt. In this paper we shall calculate 
instantons in the range of angles 7r/2 + 0{y/e) < 9h < — 0{y/e) at T = based on 
the imaginary time path integral method and their contribution to the quantum nucleation 
exponent 5q, and discuss the 9h dependence of Bq and the crossover temperature Tc from 
quantum nucleation to thermally assisted process. 



Our considerations begin by writing the transition amplitude in the imaginary time t = it 
from —Tj/2 to tj/2 based on the spin coherent state path integral formalism 

{z\exp{-TjH/h)\z) = J D[M{r,T)]exp{-SE/h), (1) 

where z is taken to be an initial easy direction of the magnetization, and Se is the Euclidean 
action which includes the Euclidean Lagrangian density L^; as 

Se= fdrf dhlE. (2) 



Since the magnitude of the total magnetization M is frozen up at sufficiently low temper- 
ature, the direction of the total magnetization becomes the only dynamical variable which 
we are interested in. In this situation the measure of the path integral i^[M(r, r)] in Eq. 
(15) is replaced by DVL as 

N 2 / + 1 

Dn = lim n ) sin 9^ d9k d<Pk , (3) 



where e = maxijk+i — Tk) and J = Mo/h'y. Here 7 = g^s/^, with g being the g-factor, /i^ 
the Bohr magneton, Mq the magnitude of the magnetization. Using the standard instanton 
method, the rate of quantum nucleation has the relation 

Tq oc exp{-SE/h), (4) 

and its Euchdean action for this problem is given by 

Se= f{iM^[i-cose{r,T)]^^^j^ + E[M{r,T)]}dTdh. (5) 
J 7 dr 

The first term in Eq. (R) is the topological Wess-Zumino term P], and the second term is the 



energy density which is composed of the magnetic anisotropy energy, the exchange energy 
and the energy given by an external magnetic field H. So, the energy taken into account 
can be written as 

E[M{r, t)] = Ea + Ee.-M-U. (6) 



We notice that from 6Se = the action (w produces the classical equation of motion |10 



whose solution is called instanton, |11| and describes the 1 © 1 dimensional dynamics in the 



Hamiltonian formulation, which consists of the canonical coordinates (p and = 1 — cos 6. 
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In this paper we study the biaxial symmetry whose magnetic anisotropy energy is given 

by 

Ea = K,{al + al) + K2al (7) 

where Ki and K2 are the parallel and transverse anisotropy constants, and ai,a2, and 
the directional cosines. The exchange interaction is expressed as [0] 

E,, = ^C[{Va,r + {Va^r + (Vas)^] (8) 

where C is an exchange constant. Applying the magnetic field in the xz-plane ||^ and ex- 
pressing «!, a2 and in terms of the fields ^(r, r) and 0(r, r) which are spherical coordinates 
of M, the total energy @ is given by 



E[e{r, r), 0(r, r)] = sin^ 6 + K2 sin^ sin^ 6 + ^C[(V^)^ + (V0)^ sin^ ^] 

-H^Mosmecos(f)- H.Mocose + Eo, (9) 

where i?o is introduced to make E{6, 0) zero at the initial orientation. As will be seen later, 
the effective mass of the magnetic particle is inversely proportional to a linear combination 
of K2 and Hx while there is no exact analog of the kinetic energy in the action (^. For 
this reason we need either the transverse anisotropy constant K2 or the magnetic field Hr^ 
transverse to the initial easy axis for quantum nucleation. Introducing the dimensionless 
constants, 

K2 = K2/2Ki, C = C/4Ki, H, = HjHo, H, = HjHo, (10) 

where Hq = 2Ki/Mo, we obtain the total energy written as 

^[^(r, r), 0(r, r)] = i sin^ 6 + K2 sin^ sin^ 6 + C[{Vef + (V0)2 sin^ 6] 

—H^sinO cos(f) — HzCosO + Eq, (11) 

where E{6, 0) = E{6, (f))/2Ki. The plane given by = is the easy plane, on which E{6, 0) 
is given by 

E{e, 0) = - sin^ e-H cos{e - Oh) + Eq, (12) 

where we used = H sin 6h and = H cos Oh since H lies in x2;-plane. Let us define 
6*0 to be the angle of metastable state generated by the anisotropy energy and the exter- 
nal magnetic field, and 9c the angle at which the barrier vanishes by the applied critical 
magnetic field ifc- Then, 6^0 is determined by [dE{6,Q)/d6]o=o^^ = 0, and 9c and He by 
[dE{9,0)/d9]()^g^^H=H, = and [d'^E{9,0)/d9\^e^^H=H, = 0. After a little calculations, the 



dimensionless critical field and the critical angle 9c are expressed as [ITS 

He = (sin2/3 9h + I cos9h\'^T''^\ (13) 
1 
2 



2| cot 9h 


1/3 


1 + 


cot 9h 


2/3 



9c = I arcsin[ ^JJ^' ]. (14) 



For, e.g., Oh = 37r/4 and vr, we have 6c = vr/4 and 0, respectively. Simple analysis for Eq. 



(|T^) shows that, since the height and the width of the barrier are proportional to the power 
of the parameter e, the value of e should be small in order to have a large tunneling rate. 
The small value of e can be achieved by tuning H to He- From now on we will consider the 
value of e to be small, that is to say, e <^ 1. From the equations which Oq and 9c satisfy, we 
obtain the equation for t]{= 9c — 9q) as ||13| 



sin(2^,) ie~^ri^)-ri cos(2^^e) (2e - r]^) = 0. (15) 
Simple calculations show that 1] is of the order of ^/e. Thus the order of magnitude of the 



second term in Eq. (|T5|) is smaller than that of the first term by ^/e and the value of rj is 



determined by the first term which leads to rj ~ \j2e/'i. However, when 9h is very close 
to 7r/2 or vr, sin(2^^c) becomes almost zero, as shown in Fig. |l], and the first term is much 
smaller than the second term in Eq. (|15]). Therefore the value of r] is obtained from the 
second term when ^^/^ ~ 7r/2 or tt, which leads to ~ \/2e for 9h — 7r/2 and r/ ~ for ~ vr 
from the detailed analysis of the equations which 9q and 9c satisfy. Since the first term in 



(P^) is dominant in the range of value 9c which satisfies tan(26'c) > O(v^), r] ^ Y2e/3 is 
valid for 7r/2 + O(y^) < 9h < tc — 0{y/e) by usnig Eq. (0). As is shown in Fig. ||, this is 
checked by performing the numerical calculation for the equations which 9o and 9c satisfy. 
From Eq. ( pT]) we obtain an approximate form of E{9, (p) as 

E{6, (j)) = K2 sin^ (f) sm\9o + 6) + sm{9o + 5) (1 - cos 0) 

+C[{V6)^ + {V<j)f sm\9o + 6)] + Ei{6), (16) 

where E{9, 0) is written as E{6, (p) by introducing a small variable 5 = 9 — 9q, and Ei{5) is 
a function of only 5 given by 

m) = \ smi29c)i35'v " 5') + ^ cos(2^^,) [^^(e - ^^') + - j]- (17) 

As previously considered for 77 in Eq. (|T^), even though the cos(26'c)-term in Eq. ( P^ ) looks 
smaller by a factor of r] which is of the order of ^e, the second term can not be neglected 



near Oh = vr/2 and vr because sin(2^c) is almost zero for these regions of Oh- In order to 
evaluate the order of magnitude of each term in the action and simplify the calculations 
for e <^ 1, it is convenient to use new scaled variables 

f = et cuor, f = r/ro, 5 = tq = V2C, ooo = 2^KjMo, (18) 

where it is noted that tq is the thickness of the domain wall in a bulk ferromagnet. Then, 
the Euchdean action (H) becomes 

SEm,f),(j){i,r)]=nje-^rlJ dfd'fite'^il-cosiOo + ^S)]^ 

+i?2sinVsin^(6'o + \/e6) + H^sm{Oo + y/e6){l - cos0) 
+ lei+^(W)2 + ^ef^sin^Oo + v^5)(V0)2 

+ \ sin(2^,)ei(35^^ " + ^ cos(2^^,)e^[5^(l - ^) + - j]}. (19) 

where the parameters a and (3 are fixed by the analysis of the order of magnitude of each 
term in the action ([T9|) subject to the situation considered. 

For 7r/2 + 0{y/e) < Oh < — 0{y/e), the critical angle 0^ and the angle of metastable 



state ^0 are in the range of angle from to 7r/2, and 77 ^ y 2e/3. In this situation, the action 
(|I^) becomes 



5'£;[(5(r, r),0(r, f)] = ^Je J dTdh{-ie^ sm{Oo + ^/e5) 



df 



+K2 sm^(f)sm\Oo + v^6) + sm{Oo + V^6) sm\^) 

+ le'^^iV6r + le^sm'iOo + v^5)(V0)^ 

+ ^sin(2^,)6i(35'^-5=^)}, (20) 

where we performed the integration by part for the first term and neglected the total time 
derivative. Since the last term in the action ( pOD makes the potential barrier in our problem, 
the order of magnitude of each term should be the same order as that of the last term. Since 
the order of mag nitude of each term is O(e(°+^)/20)^ ©(sin^^), O(sin2(0/2)), e^+^, O(e^02) 



and e^/^, respectively, the small values of (j) contribute to the path integral for e <^ 1, and its 
order of magnitude is expected to be e^/^. Choosing a = /3 = 1/2, the order of magnitude of 
the term which includes (V0)^ is while those of the other terms are e^/^. Thus, the fifth 
term in Eq. ( pO|) is higher and neglected. 

Performing the Gaussian integration over in the transition amplitude (|l|) and the 
measure (|^), the remaining integral is of the form 

j D[5(f,f)]exp(-5^^V^), (21) 

where from Eq. (^) the effective action is given by 

S]ff[5{v,f)]=hJrle^ j dfrfS^ [iM(|^)2 + i(W)^ + i sin(2^,)(v^5' - 5^)]. (22) 

Here the effective mass in the reduced dimension is given by 

sin 6c 

M=- ^ . (23) 

i/^. + 27^2 sin ^ ' 



Introducing the variables r = T^Jsm{26c)/M and r = fysin(2^c); we simplify the effective 
action (^) to be 

Sg^m f)] = nJrle^^^^ j d-d'A\{%f + + " (24) 

In case of a small magnetic particle of volume V Tq, S does not depend on the space 
f , which leads to the action 

S'E"[6{f)] = nJeW^Msm{29,) j df^-{^f + -^{V^f - f)]. (25) 

The corresponding classical trajectory is given by 

^"^'^'^ = cosh^(^)V4,)' (26) 



and its classical action is found to be 13 



Sci = (^j)li4^,5/4 (27) 





cot 6h 


|l/6 


v/i + t(i + l 


cot ^j^|2/3) 



The geometry of the problem in this paper is the thin film of thickness zt less than the 
size ro/e^/^ of the critical nucleus and its plane is perpendicular to the initial easy axis. After 
performing the integration over z-variable, the action (|2^) becomes 



(28) 



At zero temperature the classical solution of the action (^) has 0(3) symmetry in r and r. 
Therefore, 5 is a function of q where q = (x^ + + which the effective action (|28D 

becomes 







with the angular dependence 



1 + Icot^Hp/^ 



fieu) = / ' (30) 

V2\ cot - e + 1^(1 + I cot 



The corresponding classical equation of motion satisfies 

^ + (31) 

dq^ qdq 2 A ^ ^ 

This equation can be solved by the numerical approach, whose results are similar to the 
ones in Ref. 0. The maximal value of 6 is 6.85 at f = and f = 0. Using the results, we 
obtain the rate of quantum nucleation given by 

Tq oc exp(-5^/;i) =exp[-105.2Jr2/(^^)z,e=^/^]. (32) 

Since the period tj of the instanton becomes zero for the temperature much larger than 
Tc, the instanton solution becomes independent of f. In this case we have the rate of thermal 
nucleation 

f/e// 

Tt oc exp(— -— ) (33) 
where the effective height of the barrier 



t/^// = 4vr K,rl z,e f + - \f]. (34) 

From the saddle point of the functional the shape of thermal nucleus is determined by 

^,lf_fj,2?.0. (35) 

dq^ qdq 2 A ^ ^ 

The solution is found by the numerical method whose shape is similar to the one in Ref. 
0. The maximal size of the thermal nucleus is 3.9 at r = 0. Using this result, the effective 
height of the barrier becomes 

U^ff = A1.3Kirlzte. (36) 

Comparing this with Eq. (0), we obtain the approximate form of the temperature of the 
crossover from thermal to quantum regime given by 

Here we notice that the angular and the field dependence of coincides with the one of 



Ref. and with that of Ref. in the limit of K2 = 0. It is understood from the fact 
that Tc ~ hoob/ks, where the frequency of small oscillations around the minimum of the 
inverted potential is proportional to uj^e^^'^ / /(Oh) from Eq. (|18]) and the scaling relation 
between f and f . 

Let us now consider the situation with the magentic field opposite to the initial easy 
axis. For 6h = tt, we obtain 6^ = and r/ = by using Eqs. (|1^) and ([T5|). In this case, 



noting that ^0 = 0, = and e ^ 1, the Euclidean action ([T9|) can be simplified to be 

Se^, f), 0(f , f)] = hJe'^rl J dfd't[^e'+^ 6' ^ + K,e f sin^ 



+ + \e'^^6\v<pf + le\6' - ^)]. (38) 

In order to obtain the quantum tunneling under the barrier, we need the second term 
which includes the transverse component of the anisotropy energy Ea- Without K2 in Ea, 
Mz, as a quantum operator, commute with Ea and conserves. Thus, K2 is responsible for 
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quantum transition at Oh = vr. Comparing the second term with the last term which makes 
the potential barrier, only small values of cj) contribute to the path integral because e ^ 1. 
Thus we expect the order of magnitude of to be e^/^. Choosing a = /5 = 1, the order 
of magnitude of each term except the fourth term in Eq. ( pH]) is while the fourth one 
is which is higher and neglected. Following the same procedure as we have done pre- 
viously, we obtained the WKB exponent Bq{T = 0,9h = vr), the effective height of the 
barrier U'^^^ [Oh = vr) for T > Tc and the approximate form of the crossover temperature 
Tc, which were presented by Chudnovsky and Gunther f^. Recently, Ferrera and Chud- 
novsky discussed the thermally assisted quantum nucleation a.t 6h = tt and obtained the 
temperature dependence of the WKB exponent B{T, Oh = it) in the range of temperature 
< T < 2Tc by numerically solving the partial differential equation with period fiuje^/"^ /T 



where uj = 

As summarized in Table |, Bq is proportional to e^/^ for a range of angles '^^/2 + 0{^/e) < 
9h < TT — 0{^/e), while Bq oc e^/^ at 9h = tt. By using the interpolation method, we obtain 
Bq in the range of angle close to vr, and present the results in Fig. ^ for e = 0.01. In 
case of volume V <^ Bq can be minimum aX 6h = depending on the ratio 

K2/K1 of the anisotropy constants. Bq{6h = tt) is minimum for K2 = Ki while it is 
maximum in the limit of K2 <^ Ki which is uniaxial symmetry |Q. In the ferromagetic 
film which we have considered with small e, Bq has a minimum in the range of angle 
7r/2 + 0{^) < 6h < n - 0{y/e) for both K2 = Ki and K2 < Ki, as is shown in Table 
|. Also, in the limit of K2 ^ Ki the dependence of /{Oh) on the ratio K2/K1 is the same 
as that of Bq at Oh = tt. Thus, in the small e limit, the ratio does not make an effect 
on the shape of Bq, which is different from the situation for V r^. [jl3[ In other words, 
since Bq has a shape similar to Fig. ^, it always has a minimum in the range of angle 
7r/2 + 0{^/e) < Oh < — 0{^/e). The minimal value of Bq or equivalently the maximal 
value of Fq at T = can be obtained by the analysis of f{OH) and its Oh is given by 

2K 

Oh = 71- arctan[^= ^ 1^/^, (39) 

^(2^ + |P + 2-| 
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where we notice that 9h ~ 1.97(^ 112.9°) for, e.g., Ki = K2, as is shown in Fig. |^. Simple 
analysis shows that the angle which gives the minimal valus of Bq, Eq. (^), decreases as 
K2 decreases, and converges to tt — arctan(3y3)(^ 1.76 ~ 100.9°) in the limit of K2 = 0. 
1^ Since the effective height of barrier W^-^^ is proportional to e for both 7r/2 + 0{^/e) < 
9h <7^ - 0{^/e) and 9h = tt, we obtain oc e^/^ for n/2 + 0{^/e) < 9h < n - 0{y/e) and 
= IT. Since is inversely proportional to f{9H) for a wide range of angles, 
Tc becomes maximum at the angle (|39|) . For Ki ~ K2 — 10^ erg /cm?, Mq ~ 500 emu/cm^ 
and e ~ 10^^ — 10~^, the maximal value of Tc is approximately 235 — 132 mK while Tc ~ 
167 — 53 mK for 6'/^ = vr. 

In conclusion, we have studied the quantum nucleation of a thin ferromagnetic film 
placed at some angle with the magnetic field and presented the quantum nucleation exponent 
and the crossover temperature based on the spin coherent path integral formalism. If the 
experiment is to be performed, there are three control parameters for comparison with 
theory: the angle of the magnetic field 9hi the magnitude of the field H in terms of e, and 
the temperature T. The dependence of Bq and Tc on the power of e seems to be observed 
by changing 9^. Also, the angle (BUI) which gives the maximal value of Tc as well as Fg is 
expected to be observed in experiment. 
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TABLES 

TABLE L Summary of the results for the quantum nucleation in the range of angle 
7r/2 + 0{^/e) < 0H < IT - 0{^/e) and Oh = vr. Here Bq is the WKB exponent at T = and 
fien) given by Eq. @). 



Field angle 


tt/2 + O(v^) <eH<7r- O(v^) 


^H = vr i 


a = (3 


1/2 


1 


M 




(47^2)-! 


Bq/Jrlzt 


105.2 fieH)e^/' 


37.9 ^Ki/K2 




41.3 e 


23.4 e 


{kBT,)J/Ki 


0.22 e'/yfiOn) 


0.62 ^K2/Kit^'^ 
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FIGURES 

FIG. 1. Comparison of (a) sin(2(9c) with (b) cos(2(9c) in Eqs. and (||). Note that 

sm{29c) = 0, I cos(2^c)| = 1 for both 6h = vr/2 and vr. 

FIG. 2. The shape of r/(= 9c - Oq) as a function of Oh for (a) e = 0.01 and (b) e = 0.001. Note 



that 7] y/2e/3 is vahd in a wide range of angles, i.e., 7r/2 + 0{^) < 9}{ < it — 0{^/e). 

FIG. 3. The 9h dependence of the quantum nucleation exponent Bq for Ki = K2 and e = 0.01, 
where Bq = Bq/ Jr^zt- The minimal value of Bq is approximately 2.65 at 9h ~ 1.97 while 
Bq k, 3.79 at 6*/^ = vr. 



15 



0.16 



0.14 - 



0.12 - 



0.1 - 



0.08 - 



0.06 - 



<> 

<> 
<> 

<> 

<> 
<> 



(a) 
(b) 




■^-^-^-^0 C C V- 



-0 0- 




<> 



0.04 - 



0.02 - 







+ 



1.4 1.6 1.8 



2.2 2.4 2.6 2.8 



Oh 



<> 
<> 
<> 



-^-H 



<> 

<> H 

<> 
<> 



3.2 




Oh 



